Phase separation and flux quantization in the doped quantum dimer model on the 

square and triangular lattices 
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The doped two-dimensional quantum dimer model is investigated by numerical techniques on the 
square and triangular lattices, with significantly different results. On the square lattice, at small 
enough doping, there is always a phase separation between an insulating valence-bond solid and 
a uniform superfluid phase, whereas on the triangular lattice, doping leads directly to a uniform 
superfluid in a large portion of the RVB phase. Under an applied Aharonov-Bohm flux, the superfluid 
exhibits quantization in terms of half-flux quanta, consistent with Q = 2e elementary charge quanta 
in transport properties. 
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Understanding electron pairing in high temperature 
superconductors is a major challenge in strongly corre- 
lated systems. In his milestone paper, Anderson pro- 
posed a simple connection between high temperature su- 
perconductors and Mott insulators Electron pairs 
"hidden" in the strongly correlated insulating parent 
state as Valence Bond (VB) singlets lead, once fried to 
move at finite doping, to a superconducting behavior. 
A very good candidate of the insulating parent state is 
the resonating VB state (RVB), a state with only expo- 
nentially decaying correlations and no lattice symmetry 
breaking. A simple realization of RVB has been pro- 
posed by Rokhsar and Kivclson (RK) in the framework 
of an effective quantum dimer model (QDM) with only 
local processes and orthogonal dimer coverings Q ■ Even 
though the relevance of these models for the description 
of SU(2) Heiscnbcrg models is still debated, this approach 
is expected to capture the physics of systems that nat- 
urally possess singlet ground states (GS). For instance, 
specific quantum dimer models have recently been de- 
rived from a spin-orbital model describing LiNi02 [H , or 
from the trimerized kagome antiferromagnet Q . In a re- 
cent work, a family of doped QDMs (at T=0) generalizing 
the so-called RKpoint of Ref.0] has been constructed 
and investigated [5|, taking advantage of a mapping to 
classical dimer models Q that extends the mapping of 
the RK model onto a classical model at infinite temper- 
ature, with evidence of phase separation at low doping. 
However, the soluble models of Ref. [j| are 'ad hoc' con- 
structions, and this call for the investigation of similars 
issue in the context of more realistic models. In that re- 
spect, a natural minimal model to describe the motion of 
charge carriers in a sea of dimers is the two-dimensional 
quantum hard-core dimcr-gas Hamiltonian: 

H = v"£N c \c){c\-JJ2 W){c\~t £ \c")(c\ 

c (c,c') (c,c") 

(1) 

where the sum on (c) runs over all configurations of the 



Hilbert space, N c is the number of flippable plaquettes, 
the sum on (c', c) runs over all configurations |c) and |c') 
that differ by a single plaquettc dimer flip, and the sum 
on (c",c) runs over all configurations |c) and |c") that 
differ by a single hole hopping between nearest neigh- 
bors (triangular) or (diagonal) next-nearest neighbors 
(square). Throughout the energy scale is set by J = 1. 
A schematic phase diagram for the two lattices is de- 
picted in Fig[T] in the undoped case. Remarkably, these 
lattices lead to quite different insulating states. Indeed, 
an ordered plaquette phase appears on the square lattice 
immediately away from the special RK point, whereas a 
RVB liquid phase is present in the triangular lattice 
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FIG. 1: (color online) Schematic phase diagrams for the tri- 
angular and the square lattice. 

In this Letter, we investigate in details the properties 
of model (p} on the square and triangular lattices at fi- 
nite doping. Building on the differences between the two 
lattices in the undoped case, we investigate to which ex- 
tent the properties of the doped system are governed by 
the nature of the insulating parent state. This investi- 
gation is based on exact Diagonalisations and extensive 
Green's Function Monte-Carlo (GFMC) simulations Q 
essentially free of the usual finite-size limitations Q . 

Phase separation: At small t, it is expected that holes 
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experience an effective attractive potential. It is there- 
fore natural to first address the issue of phase separation 
(PS), i.e. the possibility for the system to spontaneously 
undergo a macroscopic segregation into two phases with 
different hole concentrations. We analyze the problem as 
a function of the hopping parameter t and hole concen- 
tration x = rih/N , where n% is the number of holes in the 
system and N the number of sites. In order to perform 
a Maxwell construction we define: 



s(x) 



e(x) - e(0) 



(2) 



where e(x) is the energy per site at doping x. This quan- 
tity corresponds to the slope of the line passing through 
e(0) and e(x). If the system exhibits PS, the energy will 
present a change of curvature implying s(x) to have a 
minimum at a critical doping x c . The fact that the local 
curvature of e(x) at x — is negative then implies that 
the two seperated phases will have x = x c and x = (the 
undoped insulator). In FigJ^l typical results are shown 
for both square and triangular lattices and for different 
sizes. Interestingly, PS appears in both cases, but with 
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FIG. 2: (color online) Slope of energy density (Eq.(2)) vs 
doping for different sizes, (a) Triangular lattice, (b) Square 
lattice. 

noticeable differences. While for the square lattice (lower 
panel) the critical hole concentration x c is roughly size 
independent, there is a strong size dependence for the 
triangular lattice (upper panel). This size effect can be 
traced back to the nature of the parent undoped state. 
On the square lattice, the crystalline phase (for v < 1) at 
zero doping is very robust and for increasing size, its lo- 
cal order changes only weakly. On the triangular lattice, 
it has been shown that size effects play an important role 
3, especially in the RVB liquid phase for 0.8 < v < 1. 
Periodic boundary conditions (BC) tend to stabilize the 
so-called Vl~2 x VT2 phase on small clusters, and clusters 
with more than 192 sites are necessary to significantly re- 
duce finite-size effects, in particular, as in Fig[21 close to 
the transition point with the crystalline phase. Hence the 



PS observed around x = 0.075 for the 3x6x6 cluster is 
not representative of the thermodynamic limit. 

To obtain the phase diagram in the (v, x) plane, we 
have performed a systematic size-scaling analysis at fixed 
t and for various v's depicted in FigO In agreement with 
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FIG. 3: (color online) Scaling of the critical doping x c de- 
fined by Maxwell construction with the inverse total number 
of sites, (a) Triangular lattice, (b) Square lattice. 

the previous discussion, a significant size dependence is 
only present for the triangular lattice, in which case PS 
disappears for large clusters in the RVB phase in the 
vicinity of the RK point In FigHl we report the ther- 
modynamic limit of x c for the two lattices as a function 
of v, and for different values of t. For the square lattice, 
calculations have been done from the RK point down 
to the expected phase transition between the plaq uette 
phase and the columnar phase, namely v ~ 0.6 [10j. For 
the triangular lattice, the range between the RK point 
down to the RVB--\/l2 x \/l2 transition point at v ~ 0.8 
has been covered (§]. These results clearly demonstrate 
the difference between the square and triangular lattices. 
In the first case, as soon as v 7^ 1, PS occurs for x < x c . 
Moreover, upon decreasing u, crystalline order strength- 
ens and, for fixed t, it is necessary to consider a higher 
concentration of holes to reach a stable conducting phase. 
Similarly, the bigger t, the lower x c . On the triangular 
lattice, a finite size-scaling analysis shows that no phase 
separation appears down to a critical value v ~ 0.9, well 
above the critical value v ~ 0.8 below which plaquette 
order sets in. Although numerical limitations prevent 
computations for smaller v and £, our results up to the 
3 x 12 x 12-site cluster provide clear evidence for a re- 
gion of PS inside the RVB region, between v ~ 0.8 and 
V ~ 0.9 p]]. 

Dimer ordering on the square lattice: Next we investi- 
gate how dimer order, known to exist at a; = 0, evolves 
under finite doping. Two sccnarii are a priori possible: i) 
the dimer order vanishes in the stable conducting phase 
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FIG. 4: (color online) Phase separation boundaries for the 
square and triangular lattices in the thermodynamic limit. 
The dashed lines correspond to the approximate location of 
the phase transition between plaquette and columnar phases 
[ic| for the square lattice and between plaquette and RVB 
phases Ql| for the triangular lattice. 









/ 




(a) 


r 




M k * 



FIG. 5: (color online) (a) First Brillouin zone, (b) Momen- 
tum dependence of the squared order parameter for the 196- 
site cluster (14x14). (c) Squared order parameter at the M 
point as a function of the hole concentration x and for differ- 
ent cluster sizes. The thermodynamic limit is depicted as a 
dashed line, with the corresponding error bar (see main text). 



immediately at x c ; ii) dimer order survives above x c in 
a narrow region of the conducting phase. To solve this 
problem, we have calculated the squared order parameter 
D 2 (k) in the GS |* ) defined by: 



D 2 (k) = 



1 <* |d(-*)d(fc)|*o) 



N 



(*o|*o) 



(3) 



along the path of Figj5]Ja) the first Brillouin zone of the 
square lattice, where d(k) is the Fourier transform of the 
dimer operator defined on the horizontal bonds. Note 
that this calculation has not been tried for the triangular 
lattice since no Bragg peak is present in the RVB phase, 
and the algorithm is losing efficiency for v < 0.8 [9(. In 
the pure plaquette phase on the square lattice, a Bragg 
peak develops at point fcjvf = (""j 0), the middle of the side 
of the BZ. We show in Fig[5jb) a typical result for the 
squared order parameter on the 196-site cluster for differ- 
ent values of x. Clearly, the Bragg peak disappears upon 
doping. A finite size scaling of the order parameter can 
be performed thanks to the linear behaviour at low con- 
centration. Within our data, D 2 (k = k M ) = D 2 K[ {L,x) 
behaves like a^x + in the linear region. In this 
case, one can determine rather precisely x +00 such that 
Dm (+oo, x+oo) = i.e. ir+oo is the concentration in 
the thermodynamic limit where the Bragg peak vanishes. 
By definition, x +QO = —b +00 /a +00 ~ 0.05(8). This value, 
and the linear behaviour in the thermodynamic limit, are 
displayed in FigO^c). If we compare x+oo to the corre- 
sponding x c from FigHl we can conclude that the Bragg 
peak indeed vanishes at x c ~ 0.067(5) within error bars. 
Note that numerical errors increase for larger clusters, 
and we are not able to use the same analysis for clusters 
larger than 256 sites (the results for the 18 x 18 cluster 
have not been used). Although the determination of ie+oq 



is delicate, the linear behavior of D 2 vs x is consistent 
with the physical behavior expected for the binary sys- 
tem of FigJU No dimer order is present above x c , showing 
that the system is simply "conducting "in this case, with 
D 2 (kn,j) decreasing as .ZV -1 (critical behaviour). 

Flux quantization: Finally, let us turn to a better char- 
acterization of the "conducting phase" . Since holes are 
bosonic one expects the conducting phase to be a su- 
perfluid (through Bose condensation) [l3( . However, ex- 
tra complexity results from the fact that these bosons 
move in a dimer environment. To investigate this is- 
sue, we pierce the torus by an Aharonov-Bohm flux 
of strength <f> = £0o with < £ < 1 and where 
4>q = hc/e is the magnetic flux quantum. To reduce 
finite size effects, we also consider arbitrary BC in the 
second direction (y). All this is implemented by the 
Peierls substitution, changing the hole hopping term into 
t 1 = cxp(±i2Tr£ i a/ L x ±i2irna/ L y )t, where the ± depends 
on the directions ±x±j/, while a is the lattice parameter 
and L x and L y the linear sizes of the system. Obviously, 
the whole spectrum should be periodic in £ with period 
1. We show in Fig O the spectrum of the 4x4 cluster 
on the square lattice, with 4 holes. It turns out that (i) 
the GS energy exhibits well-defined minima and (ii) is 
rigorously periodic with period £ = 1/2, which means 
that there is flux quantization in units of half the flux 
quantum (red curve). Property (i) is typical of a super- 
fluid [3] : It is the precursor of a finite barrier in the ther- 
modynamic limit. [l5j. Property (ii) was suggested quite 
some time ago in the context of a more general QDM 
by Kivelson[l3|, who also predicted that, in the cylinder 
geometry, one should be able to tunnel between the two 
branches of Fig|6j thus lifting the degeneracy at the level 
crossing. This degeneracy is not lifted in our case, neither 
in the torus geometry, nor in the cylinder geometry, due 
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FIG. 6: (color online) Energy spectrum vs (reduced) 
Aharonov-Bohm flux f, for a 4 x 4 cluster with 4 holes at 
v = 0.70 and t = 1.00. (a) Torus geometry. Arbitrary BC are 
used in the transverse direction to obtain a continuous set of 
momenta leading to a continuum (colored area) . (b) The two 
low-energy branches for the case of a cylinder (dashed line) 
and including a small bond disorder (full red line). 



to the translational symmetry, which puts the two states 
that are degenerate at £ = 1/4 into different symmetry 
sectors. However, getting rid of the translational sym- 
metry by changing the amplitude of a local dimer flip 
indeed removes the degeneracy (upper panel of FigJS]), 
leading to a detectable flux quantization in units of half 
the flux quantum in an experiment in which the flux is 
swecped. Thus, in our model, the ground-state energy 
has periodicity hc/2e, consistent with mobile elementary 
particles of charge Q — 2e in the system. Unlike what 
was recently found in a bosonic model with correlated 
hopping (ljj. these particles are not boson pairs: ^From 
the bosonic point of view, it is the statistical flux of the 
dimer background that leads to the half-flux quantiza- 
tion. If dimers are interpreted as SU(2) electron singlets, 
these singlets are the physical pairs that lead to half-flux 
quantization. This scenario is fundamentally different 
from the usual mechanism related to real space pairing 
of the charge carriers found e.g. in the extended Hubbard 
chain 17 1, in the 2-leg ladders [l8| or, more generally, in 
Luther-Emery liquids (l9j . as can be inferred from the 
exact degeneracy between £ = and £ = 0.5 for finite 
systems in the present case, to be contrasted with the 
significant finite-size effects of the other cases. 

Summary and conclusions: The numerical investiga- 
tion with Green's function Quantum Monte Carlo and ex- 
act diagonalizations of the doped two-dimensional quan- 
tum hard-core dimer model on the square and triangular 
lattices has led to a number of interesting conclusions re- 
garding hole motion in a dimer background. Phase sep- 
aration is often present at low doping, as suggested by 
earlier investigations, but our results indicate that it is 
related to the presence of valence bond order [2(| : In the 



RVB phase of the triangular lattice, PS only occurs close 
to the plaquctte phase, where short-range dimer correla- 
tions are already strong enough. Close to the RK point, 
doping the RVB phase leads directly to a superfluid phase 
as shown from its response to an Aharonov-Bohm flux. 
Moreover, we observed that the flux quantization is in 
units of half a flux quantum, consistent with the idea 
that the dimer background leads to effective particles of 
charge 2e. All these results are in qualitative agreement 
with the gauge theories of high T c superconductivity in 



strongly correlated systems [21 1 
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